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ABSTRACT 


Recent measurements of the magnetic field in the vicinity of the 
moon by Explorer 35 indicate the absence of a bow shock wave and only 
small perturbations in the lunar wake. A theoretical model has been 
developed based upon a guiding-center approximation. A favorable 
comparison of the theoretical and the experimental results is obtained. 



X. INTRODUCTION 


The solar wind represents the supersonic expansion of the 
solar corona into interplanetary space. It is comppsed primarily 
of protons and electrons with a small (less than 107.) admixture of 
doubly ionized helium and smaller amounts of heavier elements in 
various states of ionization. At the orbit of earth the solar 
wind velocity varies between 300 and 800 km/sec with a typical 
density of 5 protons /cm^. This represents a proton energy of from 
400 eV to 4 keV and a flux of approximately 3 x 10^ protons /cm^/sec. 
Imbedded in this plasma is an extended solar magnetic field of 
3 to 8 Y (1 Y’10 oersted). In general the electron and ion 
temperatures are not equal and the thermal motion of the ions aniso- 
tropic. The temperature of the ions is approximately 10 5 °K with 
T]|/T a varying' between 1.5 and 4. The electron temperatures are found 
to be 1-3 times the ion temperatures. The length scales of interest 
in the solar wind vary from a collision length of the order of 1 AU 
to the ion gyroradius of the order of 50 km and the Debye length of 
the order of 10 meters. Both the Alfven and ion thermal velocities 
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in cislunar space are approximately 50 kra/sec so that the Mach number 
of solar wind flow is of the order of 10. 

Because qf the different magnetic properties of the moon and the 
earth, the interaction of the solar wind with these two bodies is vastly 
different. The main features of the earth's case (reviewed earlled in 
this symposium*) are: (i) the formation of a magnetosphere around which 
the bulk of the solar plasma flows, (ii) the developement of an extended 
magnetic tail in the antisolar direction, and (iii) the presence of a 
detached, collisionless bow shock wave upstream from the magnetosphere. 



The adaptation of fluid flow theory to solar wind flow past the 
earth has been quite useful in describing this interaction. It has 
been especially successful in predicting the variations of the plasma 
and magnetic fields in the boundary layer known as the magnetosheath. 

It is the purpose of this paper both to discuss recent 
experimental results obtained on the solar wind interaction with the 
moon and to present a theoretical model of this interaction based 
upon kinetic theory . 



II. OBSERVATIONS OF MAGNETIC FIELD PERTURBATIONS BY LUNAR WAKE 


The experimental results were obtained from the Explorer 35 

spacecraft which was placed into a lunar orbit on July 22, 1967 

with the following oruital parameters: aposelene = 5.4 Rm (Rm = 1738 km), 

periselene = 1,5 Rm, inclination = 169° and period = 11.5 hours. 

As they orbit the earth the moon and Explorer 35 are immersed 

repetitively in interplanetary space, the geomagnetosheath and the 

geomagnetotail. Data taken while the spacecraft is in the geomagnetic 

tail have been used to set an upper limit to the magnetic moment 
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of the moon of 10 cgs units implying a surface field of less than 4Y. 

Of primary interest in this study are the periods when the moon is 
outside the earth's bow shock in the solar wind. During such times 
measurements of the magnetic field in the vicinity of the moon have 
been made during hundreds of spacecraft orbits. The characteristics 
of a bow shock (viz. a jump in field magnitude by a factor of 2-4 and 
substantial increase in fluctuations accompanying the heating of the 
plasma) have not been observed. The interplanetary magnetic field is 
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observed to be only slightly perturbed by the presence of the moon . 

In the core of the lunar wake observations characteristically show an 
increase of the magnetic field <307. of the unperturbed value. In 
the penumbral regions, around this core, both positive and negative 
perturbations of the interplanetary field strength are often observed. 
The magnitude and geometry of the perturbations vary with conditions 
in the undisturbed solar wind, verified by simultaneous measurements 
with other satellites in cis lunar space. 



The major effect of the moon on the interplanetary medium is the 
creation of a plasma cavity 7 ; the electrons and ions impacting the 
surface are absorbed by the lunar body. There is also no evidence for 
the development of a boundary layer analogous to the earth' s magneto- 
sheath. 

A sample of magnetic field measurements in the lunar wake displaying 
relatively large perturbations of the magnetic field is shown in 
Figure 1. The spacecraft trajectory is shown projected on the ecliptic 
plane in the upper left hand corner and positionally correlated with 
the data through UT annotation. As the spacecraft enters the lunar 
wake region, an oscillating pattern of - + - anomalies is detected. 
Subsequently after egression from the core region where a positive 
anomaly some 307« above the undisturbed field magnitude of approximately 
7.5Y was detected, a set of - I- - anomalies is again observed. Note 
that the latitude angle, 9, and longitude angle, in selenocentric 
solar ecliptic coordinates, indicate the absence of appreciable 
directional variations of the fields through the lunar wake region. 

Such large variations of the magnetic field in the lunar wake are 
not always observed. Figure 2 shows an example when, without the use of 
simultaneous correlative data from another spacecraft in cislunar space, 
it would have been difficult to positively identify the umbral increase. 
The time offset shown between the two satellites is due to their 
different positions in space and the structure of the interplanetary 
medium. Here it is noticed that only a positive anomaly of approxi- 
mately 20', o is readily identified in the umbral core. Other cases 
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show even less lunar effect on the field 



These two cases are representative of the perturbations 
observed in the lunar wake. It is clear from multiple passes 
through the wake that the orientation of the interplanetary 
magnetic field, 0o , and the ratio of plasma thermal energy to magnetic 
field energy, j3, are important in determining both the geometry and 
the magnitude of the magnetic field perturbations. Interested 
readers are referred to the original papers for discussion of the 
NASA-GSFC magnetic field experiment details. 



III. THEORY 


1 . Theoretical Model 

A steady-state model for the plasma and field in the vicinity 
of the moon is considered here. Let Uq denote the solar wind 
velocity in the unperturbed region upstream from the moon, Bq the 
unperturbed magnetic field, and 0 O the direction angle (the angle 
between -Uq and B Q ). The ratio of U Q to the upstream ion thermal 
velocity is of the order of 10. The surface of the moon is assumed ■ 
to absorb all impinging charged particles. 

Because the Larmor radius of charged particles is small (<3%) 
compared with the dimension of the moon, each particle may be 
approximated by a guiding- center particle 1 **. The lowest-order guiding- 
center gas has thermal motion only along the field lines^, and thus 
under this approximation, the problem is reduced to the interaction 
between a sphere and a supersonic one-dimensional guiding-center gas. 

The parallel thermal motion is assumed to be maxwellian for the guiding- 
center gas upstream, and the magnetic moment is invariant following each 
guiding- center . 

The velocity of a charged particle, V, may be expressed as the 
vector sum of the drift velocity u, the component of the guiding-center 
velocity parallel to the field v, and the Larmor gyrational velocity s_: 

I = u + v + 1 • 

To lowest order, s is related to the magnetic moment by 

p. =* % m s^/B. 



If a guiding-center distribution function, f, is defined as a function 

of v, p., and t, then f is related to the ordinary distribution 
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function F (V, x, t) by 

F ~ m f/2TTB. 

The kinetic equation can be reduced to state that f/B is conserved 
following the trajectory of a guiding center in phase space^ # 

D ( f/B) = 0. (1) 

Along each trajectory the plasma is microscopically frozen into the 
field. In the interaction region, the guiding-center distribution 
function is zero if the guiding-center trajectory is intercepted by 
the surface of the moon. Otherwise, the frozen- in condition gives 

f ( v, jj,, x, t ) = ( B/B q ) f ( v o , n, Xo, t ). (2) 

2. Analytical Solutions for Ion Flow 

In the interaction region, the guiding-center trajectory of ions 
is deflected slightly from a straight line by variations of the magnetic 
and the electric fields. Simple analytical results for the ion flow 
in the interaction region are obtained in this paper by approximating 
the ion guiding-center trajectories by straight lines. Included in 
the analysis are (i) the truncation of the distribution function due 
to absorption of particles by the surface of the moon, and (ii) the 
microscopic frozen- in condition of the guiding-center plasma. This 
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paper presents a modification of the earlier model * by including the 
frozen- in condition. 

Under the assumptions stated above, one can calculate the guiding- 
center distribution function for ions and then calculate the macroscopic 
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flow properties . The coordinate system is shown in Figure 3, the X axis 

is parallel to Ho* the Z axis is parallel to Bo x and the origin i 3 

located at the center of the moon. The scale of the coordinate system 

is normalized by using the radius of the moon as a unit length. The 

interaction region is thus intrinsically divided into four regions as 

shown in Figure 3. Let sin 2 Y ** (1 - Z 2 )/(X^ + Y^), tan X » Y/X and 
2 2 

S = U Q /(2kT||o/ni£). Then one can integrate the guiding-center 
distribution function and obtain the following results. In region 0 

n s N/N q = B/Bo (3) 

p 2 P||/P|io = B/B c (4) 

In region 1 

n = (X/2) (B/B 0 ) erfc (Yi S) (5) 

P ** n+ S (B/B q ) Yi expC-Yi 53 S 3 )//^ (6) 


where 


Yj = sin (u + X)/sin(<^o - or - X) , 


and erfc is the complement ary error function defined as 



In region 2 


erfc (x) = 1 


-t 2 

erf (x) ~ 1 - /rr*Jo e 


n = (1/2 ) (B/Bq ) erfc (YgS), (7) 

y » n + S (B/Bo)V 2 exp(-Y| S 2 )/Jn~, (8) 


where 


Y 2 sin(a - X)/sin(<i? 0 -Kv-X) . 


In region 3 


n = (1/2) (B/B 0 ) [erfc(YiS) + erfc (Y a S)] 



(9) 
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3 . Maxwell’s Equations and Numerical Solutions 

Under the assumptions that the total electric current in the 

lunar wake 3 4 * * * * * * is composed of the magnetization current, the gradient 

drift current and the curvature drift current, Maxwell's equations can 

be written in the foiloxring dimensionless form 11 , 

V . b_ = 0 (11) 

(n + iiil) V_xb_=b_xV_n •! T] p e^ x (e^ • V. ), (12) 

where £ = ^/B 0 , ej=b/b, 3= 8 Pxo/Bo 2 , and T> = T || 0 /Tjlo • 





In this frozen- in model, n and p are functions of b, x, S and % . 
Substituting the analystical solutions for n and p into Equations 
(11) and (12), one can study their solutions for n, p and as functions 
of x and four dimensionless parameters: T|, S, 0 O and 0, Following 
the method of solutions described in Reference 11, one can calculate 
the numerical solutions of Equations (11) and (12) on the XY-plane. 
Quantities obtained in a typical numerical solution are plotted as 
functions of the longitude angle, 0, at fixed radial distances in 
Figures 4 and 5. The values used for the four dimensionless parameters 


are also shown. 


IV. SUMMARY 


Experimentally the following features of the field and plasma 
In the wake of the moon are observed: (1) The magnitude of the 
magnetic field Increases In the optical shadow and a pattern of 
alternating decreases and increases is detected at times outside 
the shadot^* 6 .. (11) A plasma cavity Is produced in the core region 
of the optical shadow 7 where the field is always observed to be 
increased. (Ill) The frozen-in condition leads to positive 
correlations between the ion density and magnetic field strength, 
i*e. Jj£ « outside the optical shadow* 3 . The kinetic theory model 
discussed herein is generally capable of predicting these features. 
However, the theoretical results outside the umbral region yield 
perturbations of the magnetic field which are sometimes more than a 
factor of two smaller than the observed changes. 

In addition, the present theory also predicts that only those 
particles with high parallel thermal energy can reach the umbral region, 
while the perpendicular thermal energy of particles is directly 
proportional to field intensity. Thus the plasma is predicted to be 
highly anisotropic in the lunar wake. 

The explicit effects of electron flow and electric fields have 
been omitted. Since the flow is highly subsonic for the electrons, one 
expects their effects to be equivalent to use of a lower equivalent 
speed ratio S as they try to fill in the wake faster than the ions. 
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FIGURE CAPTIONS 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Figure 6 


Magnetic field measurements in lunar orbit obtained on 
August 7, 1967, by Explorer 35. 

Simultaneous Explorer 33-Explorer 35 magnetic field 
measurements obtained on September 10, 1967. 

Four characteristic regions of perturbed plasma flow in 
the vicinity of the moon. 

Theoretical results on the plane of symmetry for 
R/Rm - 3. 

Theoretical results bn the plane of symmetry for 

R/Rm ■ 6. 

A typical plot of the magnetic field topology on the plane 
of symmetry in the lunar wake showing the magnetic field 
magnitude perturbations and the plasma umbra. 
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